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Abstract 

The results of the HEIDELBERG-MOSCOW experiment which searches with II kg 
of enriched 7e Ge for double beta decay in the GRAN Sasso underground labo- 
ratory are presented for the full running period August 1990 - May 2003. The 
duty cycle of the experiment was ~80%, the collected statistics is 71.7 kg y. The 
background achieved in the energy region of the Q value for double beta decay is 
0.11 events/ kg y keV. The two-neutrino accompanied half-life is determined on the 
basis of more than 100 000 events. The confidence level for the neutrinoless signal 
has been improved to 4.2a. 



1 Introduction 



Since 40 years large experimental efforts have gone into the investigation of nuclear 
double beta decay which probably is the most sensitive way to look for (total) lepton 
number violation and probably the only way to decide the Dirac or Majorana nature 
of the neutrino. It has further perspectives to probe also other types of beyond 
standard model physics. This thorny way has been documented recently in some 
detail [1]. 
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The half-lives to explore lying, with the order of 10 years, in a range on 'half 
way' to that of proton decay, the two main experimental problems were to achieve a 
sufficient amount of double beta emitter material (source strength) and to reduce the 
background in such experiment to an extremely low level. The final dream behind 
all these efforts was less to see a standard-model allowed second-order effect of the 
weak interaction in the nucleus - the two-neutrino-accompanied decay mode - which 
has been observed meanwhile for about ten nuclei - but to see neutrinoless double 
beta decay, and with this a first hint of beyond standard model physics, yielding at 
the same time a solution of the absolute scale of the neutrino mass spectrum. 

In this paper we describe the HEIDELBERG-MOSCOW experiment, which was 
proposed in 1987 [2], and which started operation in the Gran Sasso Underground 
Laboratory with the first detector in 1990. We report the result of the measurements 
over the full period August 1990 until May 2003 (the experiment is still running, at 
the time of submission of this paper). The quality of the new data and of the present 
analysis, which has been improved in various respects, allowed us to significantly 
improve the investigation of the neutrinoless double beta decay process, and to 
deduce more stringent values on its parameters. 

We are giving in a detailed paper [3] (for earlier papers see [4,5]) a full description 
of the experimental procedure, the specific features and challenges of the data ac- 
quisition and background reduction. In this letter we mainly present the data, taken 
with a total statistics of 71.7 kgy for the period August 2, 1990 - May 20, 2003, and 
we present their analysis. 



2 Experimental Parameters 

The HEIDELBERG-MOSCOW experiment is, with five enriched (to (86% 88%)) 
high-purity p-type Germanium detectors, of in total 10.96 kg of active volume, us- 
ing the largest source strength of all double beta experiments at present, and has 
reached a record low level of background. It is since ten years now the most sensitive 
double beta decay experiment worldwide. The detectors were the first high-purity Ge 
detectors ever produced. The degree of enrichment has been checked by investiga- 
tion of tiny pieces of Ge after crystal production using the Heidelberg MP-Tandem 
accelerator as a mass spectrometer. Since 2001 the experiment has been operated 
only by the Heidelberg group, which also performed the analysis of the experiment 
from its very beginning. 

We start by listing up some of the most essential features of the experiment. 



(1) Since the sensitivity for the 0v(3(3 half-life is Tffo ~ a x ey (and ~ 
(m u )), with a denoting the degree of enrichment, e the efficiency of the detector 
for detection of a double beta event, M the detector (source) mass, AE the 
energy resolution, B the background and t the measuring time, the sensitivity 
of our 11kg of enriched 76 Ge experiment corresponds to that of an at least 
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1.2 ton natural Ge experiment. After enrichment - the other most important 
parameters of a ft (3 experiment are: energy resolution, background and source 
strength. 

(2) The high energy resolution of the Ge detectors of 0.2% or better, assures that 
there is no background for a 0iy(3(3 line from the two-neutrino double beta 
decay in this experiment (5.5 x 10~ 9 events expected in the energy range 
2035-^2039.1 keV), in contrast to most other present experimental approaches, 
where limited energy resolution is a severe drawback. 

(3) The efficiency of Ge detectors for detection of 0vf3f3 decay events is close to 
100% (95%, see [6] (a)). 

(4) The source strength in this experiment of 11kg is the largest source strength 
ever operated in a double beta decay experiment. 

(5) The background reached in this experiment, is 0. 113±0. 007 events/kg ykeV 
(in the period 1995-2003) in the 0vf3(3 decay region (around Q/3/3). This is the 
lowest limit ever obtained in such type of experiment. 

(6) The statistics collected in this experiment during 13 years of stable running is 
the largest ever collected in a double beta decay experiment. The experiment 
took data during ~ 80% of its installation time. 

(7) The Q value for neutrinoless double beta decay has been determined recently 
with high precision [7,8]. 



3 Background of the experiment 

The background of the experiment consists of 

(1) primordial activities of the natural decay chains from 238 U, 232 Th, and 40 if; 

(2) anthropogenic radio nuclides, like 137 Cs, lu Cs, 125 56, 207 Bi; 

(3) cosmogenic isotopes, produced by activation due to cosmic rays during pro- 
duction and transport. 

(4) the bremsstrahlungs spectrum of 210 Bi (daughter of 210 Pb); 

(5) elastic and inelastic neutron scattering, and 

(6) direct muon-induced events. 

The detectors, except detector No. 4, are operated in a common Pb shielding of 30 
cm, which consists of an inner shielding of 10 cm radiopure LC2-grade Pb followed 
by 20 cm of Boliden Pb. The whole setup is placed in an air-tight steel box and 
flushed with radiopure nitrogen in order to suppress the 222 Rn contamination of 



3 



the air. The shielding has been improved in the course of the measurement. The 
steel box is since 1994 centered inside a 10-cm boron-loaded polyethylene shielding 
to decrease the neutron flux from outside. An active anticoincidence shielding is 
placed on the top of the setup since 1995 to reduce the effect of muons. Detector 
No. 4 is installed in a separate setup, which has an inner shielding of 27.5 cm 
electrolytical Cu, 20 cm lead, and boron-loaded polyethylene shielding below the 
steel box, but no muon shielding. 

The setup has been kept hermetically closed since installation of detector 5 in Febru- 
ary 1995. Since then no radioactive contaminations of the inner of the experimental 
setup by air and dust from the tunnel could occur. 

A detailed description of the remaining background, and corresponding Monte Carlo 
simulations with GEANT4 are given in [9]. Important result is, that there are no 
background 7-lines to be expected at the position of an expected 0v(3(3 line, ac- 
cording to this Monte Carlo analysis of radioactive impurities in the experimental 
setup [9] and according to the compilations in [10]. Some 7-peaks in the simulated 
spectrum near are emitted from the isotope 214 Bi. So from Bi ( 238 [/-decay 
chain) in the vicinity of the Q-value of the double beta decay of = 2039 keV, 
according to the Table of Isotopes [10], weak lines should be expected at 2010.7, 
2016.7, 2021.8 and 2052.9 keV (see [11]). It should be especially noted, that also 
neutron capture reactions 74 Ge(n, 7) 75 Ge and 76 Ge(n, 7) 77 Ge and subsequent ra- 
dioactive decay have been simulated by GEANT4 in [9]. The simulation yields, in 
the range of the spectrum 1990 - 2110 keV, for the total contribution from decay of 
77 Ge 0.15 counts. So in particular, the 2037.8 keV transition in the 7-decay following 
(3~ decay of 77 Ge [10] is not expected to be seen in the measured spectra. It would 
always be accompanied by a 10 times stronger line at 2000.2 keV [10], which is not 
seen in the spectrum (Figs. 3,4). 



4 Control of stability of the experiment and data taking 

To control the stability of the experiment, a calibration with a 22S Th and 
a 152 Eu+ 228 Th source, has been done weekly. High voltage of the detectors, tem- 
perature in the detector cave and the computer room, the nitrogen flow flushing 
the detector boxes to remove radon, the muon anticoincidence signal, leakage cur- 
rent of the detectors, overall and individual trigger rates are monitored daily. The 
energy spectrum is taken in parallel in 8192 channels in the range from threshold 
up to about 3 MeV, and in a spectrum up to about 8 MeV. Data taking since 
November 1995 is done by a CAMAC system, and CETIA processor in event by 
event mode. To read out the used 250 MHz flash ADCs of type Analog Devices 9038 
JE (in DL515 modules), which allow digital measurement of pulse shapes for the 
four largest detectors (for later off-line analysis) , a data acquisition system on VME 
basis has been developed [6](c). The resolution of the FADC's is 8 bit, and thus not 
sufficient for a measurement of energies. The energy signals for high and low-energy 
spectra are recorded with 13 bit ADC's developed at MPI Heidelberg. The acquisi- 
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tion system used until 1995 (VAX/VMS) did not allow to measure the pulse shapes 
of all detectors because of the too low data transfer rates. For the period 1990-1995 
we prepared the final sum spectrum taken with the former electronics shortly after 
this measuring period, setting the proper conditions for data acceptance [6](a,b,c). 
Since 1995, in total we have taken with the since then new electronics, 2142 runs 
(10 513 data sets for the five detectors) (without calibration measurements), the 
average length of which was about one day. From these raw data runs and data sets 
only those are considered for further analysis which fulfill the following conditions: 

(1) no coincidence with another Ge detector; 

(2) no coincidence with plastic scintillator (muon shield); 

(3) no deviation from average count rate of each detector more 
than ±5a; 

(4) only pulses with ratios of the energy determined in the spectroscopy branch, 
and the area under the pulse detected in the timing branch (Eol values [12]), 
within a ±3cr range around the mean value for each detector are accepted; 

(5) we ignore for each detector the first 200 days of operation, corresponding to 
about three half-lifes of 56 Co (1^=77.27 days), to allow for decay of short- 
lived radioactive impurities; 

(6) ADC's or other electronics units were working properly (corrupted data are 
excluded). 

From the totally registered 951 044 events in the spectrum measured since 1995, 
there remain finally 786 652 events (9570 data sets). From these we find in the 
range 2000 to 2060 keV around - 562 events. 



5 Q Value 

The expectation for a 0^(3(3 signal would be a sharp line at the Q value of the process. 
Earlier measurements gave Q /3/3 =(2040.71 ± 0.52) keV, (2038.56 ± 0.32) keV and 
(2038.668 ± 2.142) keV [8]. The recent precision measurement of [7] gave a value of 
(2039.006 ± 0.050) keV. 



6 Development of the experiment 

Table 1 shows the collected statistics of the experimental setup and of the back- 
ground number in the different data acquisition periods for the five enriched detec- 
tors of the experiment for the total measuring time August 1990 to May 2003. 
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Table 1 

Collected statistics and background numbers in the different data acquisition peri- 
ods for the enriched detectors of the HEIDELBERG-MOSCOW experiment for the 
period 1990 - 2003. *) background determined by averaging counting rate in interval 
2000-2100 keV (counts/keV y kg); background determined from 11.1995 till 05.2003 
at Q^a from fit of spectrum in range 2000-2060 keV is (0.113±0.007) events/ keV y 
kg (all numbers without PSA). 



Detec- 


Life Time 


Date 


Backg- 


PSA 


Life Time 


Date 


Backg- 


PSA 


tor 


accepted 


Start End 


round*) 




accepted 


Start End 


round*) 




No. 


(days) 








(days) 








1 


930.9 


8/90-8/95 


0.31 


no 


2090.61 


11/95-5/03 


0.20 


no 


2 


997.2 


9/91-8/95 


0.21 


no 


1894.11 


11/95-5/03 


0.11 


yes 


3 


753.1 


9/92-8/95 


0.20 


no 


2079.46 


11/95-5/03 


0.17 


yes 


4 


61.0 


1/95-8/95 


0.43 


no 


1384.69 


11/95-5/03 


0.21 


yes 


5 




12/94-8/95 


0.23 


no 


2076.34 


11/95-5/03 


0.17 


yes 




Over period 1990 - 1995 


Over period 1995-2003 




Accepted life time = 


15.05 kgy 


Accepted life time = 


56.655 kgy 



7 Reliability of data acquisition and data 



We have carefully checked all data before starting the analysis. For details we refer 
to [3]. We checked the proper operation of the thresholds in all ADC's. Fig. 1 shows 
the threshold ranges for detectors 4 and 5, from the list mode data acquisition of the 
measurement, for the data used in the analysis. The effect discussed in [13] (their 
Fig. 4) is not present in our analysis. Also the other effects discussed in [13] are 
not present in the analysis. A detailed analysis shows that these effects arise from 
including corrupt data into the analysis (see also [3]). 

The stability of the experiment over the years and of the energy resolution and 
calibration is demonstrated in the following figures for the two 228 Th lines used for 
calibration in the range of Qpp: Fig. 1 (right) shows the sum of all weekly calibration 
spectra (^360 spectra for each detector) for the 2614.5 keV Th line. Fig. 2(left) shows 
the resolution and energy position after summing the spectra of the five detectors 
(~ 1800 spectra). The integrated resolution over 8 years of measurement is found 
to be 3.27 keV, i.e. better than in our earlier analysis of the data until May 2000 
[14,15]. This is a consequence also of the refined summing procedure we used for 
the individual 9 570 data sets. 

Fig. 2 (right) shows as another example the sum of all weekly calibration spectra 
for the 2103.5 keV line for each detector. The agreement of the energy positions for 
the different detectors is seen to be very good. 
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Fig. 1. Left: Threshold ranges for detectors 4 and 5 in the spectrum measured during 
1995 - 2003. Right: Sum of the (~360) weekly calibration spectra for the 2614.5 keV 
line from 22S Th over the measuring time 1995 - 2003. 

>1Q 6 x 10" 




Fig. 2. Left: Sum of the sum of the weekly calibration spectra of all five detectors 
shown in Fig. 1 (right) (~1800 spectra). Right: Sum of calibration spectra for each 
detector (absolute intensities) for the 2103.5 keV line for the period 1995 - 2003 
(from top to bottom detectors: 2,3,5,1,4). 

We also checked the arrival time distribution of the events, and proved that this 
is not affected by any technical operation during the measurement, such as the 
calibration procedure (introducing the Th source into the detector chamber through 
a thin tube) and simultaneously refilling of liquid nitrogen to the detectors etc. No 
hint is seen for introduction of any radioactive material from the exterior. The 
measured distribution corresponds to a uniform distribution in time on a 99% c.l. 
[3]. This has been proven by a Kolmogorov-Smirnov goodness-of-fit test. 



8 Data and Analysis 



Fig. 3 shows the total sum spectrum measured over the full energy range of all five 
detectors for the period August 1990 to May 2003. All identified lines are indicated 
with their source of origin (for details see [9],[6](e,f)). 
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Fig. 3. The total sum spectrum measured over the full energy range with all five 
detectors (in total 10.96 kg enriched in 76 Ge to 86%) - for the period August 1990 
to May 2003. 

In the measured spectra (Figs. 4) we see in the range around Qpp the 2lA Bi lines 
at 2010.7, 2016.7, 2021.8, 2052.9 keV, a line at Qpp and a candidate of a line at ~ 
2030 keV (the latter could originate, as mentioned in [11,16], from electron conver- 
sion of the 2118 keV 7-line from 114 Bi. On the other hand, as noted in [16], its energy 
differs from Qpp just by the K-shell X-ray energies of Ge(Se) of 9.2(10.5) keV). We 
show the energy range 2000-2060 keV, since here a simultaneous fit can be made 
to safely attributed lines (except the 2030 keV line), in contrast to the range we 
used earlier (2000-2100 keV). Non-integer numbers in the sum spectra are simply a 
binning effect. While a peak-search procedure as shown in Fig. 4 (bottom right) was 
needed earlier, to project out the lines from the spectrum [14,15], they can now, with 
more statistics and improved analysis, be directly clearly seen in the spectra. The 
large improvement of the present analysis compared to our paper from 2001 [14], is 
clearly seen from Fig. 4 (top left) showing the new analysis of the data 1990-2000, as 
performed here - to be compared to the corresponding figure in [14,15]. One reason 
is, that the stricter conditions for accepting data into the analysis, in particular of 
conditions (4) and (6) in section 4, cleaned the spectrum considerably. The spec- 
trum in Fig. 4 (top left) now corresponds to 50.57 kg y to be compared to 54.98 kgy 
in [14], for the same measuring period. The second reason is a better energy cali- 
bration of the individual runs. The third reason is the refined summing procedure 
of the individual data sets mentioned above and the correspondingly better energy 
resolution of the final spectrum. (For more details see [3]). The signal strength seen 
in the individual detectors in the period 1995-2003 is shown in Fig. 5. When adding 
the spectra of Fig. 5 for two subsets of detectors, detectors 1+2+4, and 3+5, both 
resulting spectra already clearly indicate, in contradiction to a claim of [13], the sig- 
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Fig. 4. The total sum spectrum of all five detectors (in total 10.96 kg enriched in 
Ge), in the range 2000 - 2060 keV and its fit, for the periods: Top: left - August 
1990 to May 2000 (50.57kgy); right - August 1990 to May 2003 (71.7kgy). Bottom: 
left - November 1995 to May 2003 (56.66 kgy); right - scan for lines in the spectrum 
shown on the left, with the MLM method (see text). The Bi lines at 2010.7, 2016.7, 
2021.8 and 2052.9 keV are seen, and in addition a signal at ~ 2039 keV. 

nal at Q^g (see Fig. 6 and Table 2). The time distribution of the events throughout 
the measuring time and the distribution among the detectors corresponds to the 
expectation for a constant rate, and to the masses of the detectors (see Fig. 7). 

The spectra have been analyzed by different methods: Least Squares Method, Max- 
imum Likelihood Method (MLM) and Feldman-Cousins Method. The analysis is 
performed without subtraction of any background. We always process background- 
plus-signal data since the difference between two Poissonian variables does not pro- 
duce a Poissonian distribution [17]. This point is sometimes overlooked. So, e.g., in 
[18] a formula is developed making use of such subtraction and as a consequence 
the analysis given in [18] provides overestimated standard errors. 

We have performed first a simultaneous fit of the range 2000 - 2060 keV of the 
measured spectra by the nonlinear least squares method, using the Levenberg- 
Marquardt algorithm [19]. It is applicable in any statistics [20] under the following 
conditions: 1) relative errors asymptotic to zero, 2) ratio of signal to background 
asymptotic constant. It does not require exact knowledge of the probability density 
function of the data. 
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Fig. 5. The spectra of the individual five detectors in the range 2000-2060 keV for 
the period November 1995 to May 2003. The sum corresponds to the spectrum 
shown in Fig. 4, bottom left. 
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Fig. 6. Spectra measured with detectors 1+2+4 (left part), detectors 3+5 (right 
part) over the period 1995-2003. The sum of the spectra is equal to that shown in 
Fig. 4 (bottom left). 
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Fig. 7. The measured total counting rates in the period November 1995 to May 2003 
in the energy range 2036.5-^2041.0 keV as function of time (la statistical errors are 
given) . 

We fitted the spectra using n Gaussians (n is equal to the number of lines, which 
we want to fit) G(Ei, Ej,aj) and using different background models B[Ej): simu- 
lated background (linear with fixed slope) (see [9]), linear, and constant. The fitting 
function F(Ei) thus is a sum of Gaussians and background. 

Error estimate: the Levenberg-Marquardt method [19] is one of the most tested min- 
imization algorithms, finding fits most directly and efficiently. It also is reasonably 
insensitive to the starting values of the parameters. This method has advantages over 
others of providing an estimate of the full error matrix. The MATLAB Statistical 
Toolbox [21] provides functions for confidence interval estimation for all parameters. 
We tested the confidence intervals calculated by the MATLAB statistical functions 
with numerical simulations. As done earlier for other statistical methods [15], we 
have simulated samples of 100 000 spectra each with Poisson-distributed background 
and a Gaussian-shaped (Poisson-distributed) line of given intensity, and have inves- 
tigated, in how many cases we find in the analysis the known intensities inside the 
given confidence range (for details see [3]). We find the confidence levels calculated 
by this method to be correct within ~ 0.3 a. The second standard method, we have 
used to analyze the measured spectra, is the Maximum Likelihood Method (MLM) 
using the root package from CERN [22], which exploits the program MINUIT for 
error calculation. For this purpose the program provided there had to be extended 
for application to non- integer numbers (for details see [23]). We performed also in 
this case a test of the given confidence intervals given by this program, by numerical 
simulations. Taking into account the results, the application of the MLM yields re- 
sults for the confidence levels consistent with the Least Squares Method. The third 
method used for analysis is the Feldman-Cousins-Method (FCM) described in detail 
in [24]. 



9 Results 

Fig. 4 shows together with the measured spectra in the range around (2000 
- 2060 keV), the fits using the first of the methods described. In these fits a linear 
decreasing shape of the background as function of energy was chosen, corresponding 
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to the complete simulation of the background performed in [9] by GEANT4. In 
these fits the peak positions, widths and intensities are determined simultaneously, 
and also the absolute level of the background. E.g. in Fig. 4(middle), the fitted 
background corresponds to (55.94 ± 3.92) kgy if extrapolated from the background 
simulated in [9] for the measurement with 49.59 kgy of statistics. This is almost 
exactly the statistical significance of the present experiment (56.66 kg y) and thus a 
nice proof of consistency. 

For the line near Q«g we obtain for the measuring periods 1990-2003 (1995-2003) the 
following result: energy (in keV) 2038.07±0.44 (2038.44±0.45), intensity 28.75±6.86 
(23.04±5.66) events. This result is obtained assuming the linearly decreasing back- 
ground as described. Assuming a constant background in the range 2000 - 2060 keV 
or keeping also the slope of a linearly varying background as a free parameter, yields 
very similar results. For the analysis of the other lines we refer to [3]. For the period 
1990-^2000 (see Fig. 4) we obtain an intensity of 19.36±6.22 events, consistent with 
the 14.8 events reported in [14,15] for 46.5 kgy. The signal at ~2039keV in the full 
spectrum thus reaches a 4.2 a confidence level for the period 1990-2003, and of 4.1 
o~ for the period 1995-2003. We have given a detailed comparison of the spectrum 
measured in this experiment with other Ge experiments in [16]. It is found that the 
most sensitive experiment with natural Ge detectors [25], and the first experiment 
using enriched (not yet high-purity) 76 Ge detectors [26] find essentially the same 
background lines ( 214 Bi etc.), but no indication for the line near Qpp. 

The fits obtained with the Maximum Likelihood Method agree with the results, 
given above. The estimations of the Feldman-Cousins Method yield confidence levels 
for the line at Qga at a ~ 4cr level. The uncertainty in the determination of the 
background here is larger, than in a consistent fit of the full range of interest by the 
other two methods. 



10 Time structure of events 



We developed for the HEIDELBERG-MOSCOW experiment some additional tool 
of independent verification, whose present status will be presented here briefly. The 
method is to exploit the time structure of the events and to select (5(5 events by their 
pulse shape. Applying a method similar to that described earlier [15], accepting only 
events as single site events (SSE), classified by one of the neuronal net methods [12], 
and in most cases simultaneously also by the second and third method [12] men- 
tioned in [15], as SSE (for details see [23]), we find a 3.3cr signal for the line near 
(12. 36±3. 72 events), and we find the other lines, at 2010.7, 2016.7 and 2021.8 keV 
to be consistent in intensity with normal 7-lines i.e. considerably more reduced (see 
[3]). When calibrating the pulse shape analysis method with the 1592 keV double 
escape line from 22S Th source - assuming that the 0v(5(5 events behave like the 
single site events in a double escape line - it can be concluded that the line at Qpp 
is consisting of SSE (for details see [3]). 
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1800 1850 1900 1950 2000 2050 2100 2150 2200 2250 
Energy, keV 




1800 1850 1900 1950 2000 2050 2100 2150 2200 2250 
Energy, keV 

Fig. 8. The pulse-shape selected (above) and the full (below) spectrum measured 
with detectors 2,3,4,5 from 1995-2003, in the energy range 1800-2250 keV. 

That a neuronal net method selection of the specific pulse shapes created by double 
beta decay can be done with much higher sensitivity, can be seen from Fig. 8. Here a 
defined subclass of the shapes selected above is selected, corresponding to events in 
part of the inner volume of the detectors. Except a line which sticks out sharply near 
Q/33, all other lines are very strongly suppressed. The probability to find ~7 events in 
two neighboring channels from background fluctuations is calculated to be 0.013%. 
Thus, we see a line near Qng at a 3.8<r level. This method also fulfills the criterium 
to select properly the continuous 2uf3f3 spectrum (see Fig. 9), which is not done 
properly by the earlier used PSA methods (for details see [3,23]). We consider this 
as an additional proof that we have observed neutrinoless double beta decay. For a 
detailed discusion see [23]. 




2010 2020 2030 2040 2050 

Energy, keV 



1000 1500 2000 2500 3000 

Energy ,keV 



Fig. 9. Left: The pulse shape selected spectrum in the range 2000-2060 keV, taken 
with detectors 2,3,4,5, in the period 1995-2003 (compare also to Fig. 4 bottom, left). 
Right: the pulse-shape selected spectrum, selected with the same method as in used 
in Fig. 8 (top) and in this figure, left part, and for the same measuring period, in 
the energy range of (100-^3000) keV. The solid curve corresponds to the shape of a 
2vf3f3 spectrum. 

The energy of the line seen in Figs. 8,9 is understood to be sligtly shifted to lower 
energy by ballistic effects (see [23]). 
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The 2039 keV line as a single site events signal cannot be the double escape line of 
a 7-line whose full energy peak would be expected at 3061 keV where no indication 
of a line is found in the spectrum measured up to 8 MeV (see [3]). Summarizing, the 
investigation of /3/3-like events is proving a line near consisting of 0v(3(3 events. 

Table 2 

Half-life for the neutrinoless decay mode and the deduced effective neutrino mass 
from the HEIDELBERG-MOSCOW experiment (the nuclear matrix element of [27] 
is used, see text). Shown are in addition to various accumulated total measuring 
times also the results for four non- overlapping data sets: the time periods 11.1995- 
09.1999 and 09.1999^-05.2003 for all detectors, and the time period 1995-^2003 for 
two sets of detectors: 1+2+4, and 3+5. *) denotes best value. 



Significan- 
ce [kgy] 


Detectors 


t?? 2 M 

(3a range) 


(m) [cV] 
(3a range) 


Conf. 
level (a) 


Period 1990 -=- 2003 


71.7 


1,2,3,4,5 


(0.69 - 4.18) x 10 25 
1.19 x 10 25 * 


(0.24 - 0.58) 
0.44* 


4.2 


Period 1995 -=- 2003 


56.66 


1,2,3,4,5 


(0.67 - 4.45) x 10 25 
1.17 x 10 25 * 


(0.23 - 0.59) 
0.45* 


4.1 


51.39 


2,3,4,5 


(0.68 - 7.3) x 10 25 
1.25 x 10 25 * 


(0.18 - 0.58) 
0.43* 


3.6 


42.69 


2,3,5 


(0.88 - 4.84) x 10 25 
(2a range) 
1.5 x 10 25 * 


(0.22 - 0.51) 
(2a range) 
0.39* 


2.9 


51.39 (SSE) 


2,3,4,5 


(1.04 - 20.38) x 10 25 
1.98 x 10 25 * 


(0.11 - 0.47) 
0.34* 


3.3 


28.21 


1,2,4 


(0.67 - 6.56) x 10 25 
(2a range) 
1.22 x 10 25 * 


(0.19 - 0.59) 
(2a range) 
0.44* 


2.5 


28.35 


3,5 


(0.59 - 4.29) x 10 25 
(2a range) 
1.03 x 10 25 * 


(0.23 - 0.63) 
(2a range) 
0.48* 


2.6 


Period 1990 2000 


50.57 


1,2,3,4,5 


(0.63 - 3.48) x 10 25 
1.24 x 10 25 * 


(0.08 - 0.61) 
0.43* 


3.1 


Period 1995 4- 09.1999 


26.59 


1,2,3,4,5 


(0.43 - 12.28) x 10 25 
0.84 x 10 25 * 


(0.14 - 0.73) 
0.53* 


3.2 


Period 09.1999 4 05.2003 


30.0 


1,2,3,4,5 


(0.60 - 8.4) x 10 25 
1.12 x 10 25 * 


(0.17-0.63) 
0.46* 


3.5 
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11 Half-Life of Neutrinoless Double Beta Decay of 76 Ge 



On the basis of the findings presented above we translate the observed number of 
events into the half-life for neutrinoless double beta decay. In Table 2 we give the 
half-lives deduced from the full data sets taken in 1995-2003 and in 1990-2003 keV 
and of some partial data sets, including also the values deduced from the single site 
spectrum. In the latter case we do conservatively not apply, a calibration factor. 
Also given are the effective neutrino masses. The result obtained is consistent with 
the results we reported in [14,15], and also with the limits given earlier in [5]. 
Concluding we confirm, with 4.2cj (99.9973% c.l.) probability, our claim from 2001 
[14,15] of first evidence for the neutrinoless double beta decay mode. 



12 Consequences 

Consequence of the above result is, that lepton number is not conserved. Further 
the neutrino is a Majorana particle see, e.g. [28]. Both of these conclusions are in- 
dependent of any discussion of nuclear matrix elements. The matrix element enters 
when we derive a value for the effective neutrino mass - making the most natural 
assumption that the 0u/3f3 decay amplitude is dominated by exchange of a mas- 
sive Majorana neutrino. The half-life for the neutrinoless decay mode is under this 
assumption given by [27] 

[t°; 2 (o+ - o+r 1 = C mm ^ + C m {tf + C AA (A) 2 + C mv ( V )^ + C mX (X)^ + 
C vX (r])(X), (m) = \m ( eJ\ + e i ^\m i eJ\ + e i **\mg)\ , 

where mil = \rnil | exp (i<fti) (i = 1,2,3) are the contributions to the effective 
mass (m) from individual mass eigenstates, with denoting relative Majorana 
phases connected with CP violation, and C mm ,C m , ... denote nuclear matrix el- 
ements squared. Ignoring contributions from right-handed weak currents on the 
right-hand side of the above equation, only the first term remains. 

Using the nuclear matrix element from [27], we conclude from the measured half- life 
given in Table 2 the effective mass (m) to be (m) = (0.2 4- 0.6) eV (99.73% c.l.). The 
matrix element given by [27] was the prediction closest to the later measured 2v(3(3 
decay half-life of (1.74^q']J) x 10 21 y [9,4]. It underestimates the 2v matrix elements 
by 32% and thus these calculations will also underestimate (to a smaller extent) the 
matrix element for 0vf3f3 decay, and consequently correspondingly overestimate the 
(effective) neutrino mass. Allowing conservatively for an uncertainty of the nuclear 
matrix element of ± 50%, the range for the effective mass may widen to (m) = (0.1 
- 0.9) eV (99.73% c.l.). 
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With the value deduced for the effective neutrino mass, the HEIDELBERG-MOSCOW 
experiment excludes several of the neutrino mass scenarios allowed from present neu- 
trino oscillation experiments (see Fig. 1, in [30]) - allowing only for degenerate mass 
scenarios [42]. 

Assuming other mechanisms to dominate the 0vf3f3 decay amplitude, which have 
been studied extensively in recent years, the result allows to set stringent limits on 
parameters of SUSY models, leptoquarks, compositeness, masses of heavy neutri- 
nos, of the right-handed W boson, and possible violation of Lorentz invariance and 
equivalence principle in the neutrino sector. For a discussion and for references we 
refer to [1,29]. 



13 Conclusion - Perspectives 

Concluding, evidence for a signal at Qm on a confidence level of 4.2a has been ob- 
served confirming our earlier claim. On the basis of our pulse shape analysis this can 
be interpreted as evidence for neutrinoless double beta decay of 76 Ge, and thus for 
total lepton number nonconservation, and for a non-vanishing Majorana neutrino 
mass. Recent information from many independent sides seems to condense now to a 
nonvanishing neutrino mass of the order of the value found by the HEIDELBERG- 
MOSCOW experiment. This is the case for the results from CMB, LSS, neutrino 
oscillations, particle theory and cosmology (for a detailed discussion see [3]). To 
mention a few examples: Neutrino oscillations require in the case of degenerate neu- 
trinos common mass eigenvalues of m > 0.04 eV. An analysis of CMB, large scale 
structure and X-ray from clusters of galaxies yields a 'preferred' value for m v of 
0.6 eV [31]. WMAP yields < 1.0 eV [32], SDSS yields < 1.7 eV [33]. 

Theoretical papers require degenerate neutrinos with m > 0.1 and 0.2 eV [35,36], 
and the recent alternative cosmological concordance model requires relic neutrinos 
with mass of order of eV [37]. The Z-burst scenario for ultra-high energy cosmic rays 
requires 0.4 eV [38], and also non-standard model (g-2) has been connected 

with degenerate neutrino masses >0.2eV [39]. It has been discussed that the Ma- 
jorana nature of the neutrino may tell us that spacetime does realize a construct 
that is central to construction of supersymmetric theories [40] . 

Looking into the future, there is some hope from the numbers given in Table 2, that 
the future tritium experiment KATRIN [41] may see a positive signal, if tritium 
decay can solve this problem at all (see [34]). 

From future double beta projects to improve the present accuracy of the effective 
neutrino mass one has to require that they should be able to differentiate between 
a (3 and a 7 signal, or that the tracks of the emitted electrons should be measured. 
At the same time, as is visible from the present information, the energy resolution 
should be at least in the order of that of Ge semiconductor detectors. If one wants to 
get independent additional insight into the neutrinoless double beta decay process, 
one would probably wish to see the effect in another isotope, which would then 
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simultaneously give additional information also on the nuclear matrix elements. In 
view of the above remarks, future efforts to obtain deeper information on the pro- 
cess of neutrinoless double beta decay, would require a new experimental approach, 
different from all, what is at present persued. 
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